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Abstract 

Solar, atmospheric and reactor neutrino experiments established that 
neutrinos are massive. It is quite natural then to consider neutrinos as 
candidate particles for explaining the dark matter in halos around galax- 
ies. We study the gravitational clustering of these neutrinos within a 
model of a massive core and a surrounding spherical neutrino halo. The 
neutrinos form a degenerate Fermi gas and a loaded polytropic equation 
is established. We solve the equation and we obtain the neutrino density 
in a galaxy, the size of the galaxy and the galactic rotational curves. The 
available data favor a neutrino with a mass around lOeV. The consequent 
cosmological implications are examined. 



Solar, atmospheric, reactor and accelerator neutrinos have provided com- 
pelling evidence for the existence of neutrino oscillations, implying non-zero 
neutrino masses [T]. Solar neutrino oscillations indicate a mass squared differ- 
ence Ato|q^ = 6 X 10~^eV^ while atmospheric neutrino oscillations suggest 
Am^rpj^j = 2 X 10~^eV^. Further, attributing the LSND anomaly to neutrino 
oscillations implies the existence of a fourth neutrino, a sterile one, with a mass 
above 1 eV [2]. Cosmology offers further insights on the neutrino mass scales[3]- 
[B]. It is inevitable that neutrinos have participated in gravitational clustering 
around massive galaxies, constituting part of the dark matter. Our purpose 
is to investigate the neutrino clustering effect, obtaining information on the 
cosmological presence of massive neutrinos. 

The universe today exhibits structure on many scales. Galaxies range in 
mass, with most bright galaxies having masses of 10^° — 10^^ Mq. The sizes 
and distributions of present-day galaxies reflect the spectrum of initial density 
fluctuations, seeded by random motion or cosmic strings |7]. We expect that 
each galaxy attracts the cosmological neutrinos in its neighbourhood, creating 
a neutrino halo of an average mass Mi, sa 3 x 10^" Mq. 

Our idealized galaxy consists of a spherical massive core (of a mass = 
(tMq and a radius of few kpc) surrounded by a spherical neutrino halo. 
Hydrostatic equilibrium prevails 
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where N{r) is the neutrino density, m is the neutrino mass and P{r) is the 
pressure of the neutrino gas. Considering that the neutrinos form a degenerate 
Fermi gas we obtain a polytrope rquation 
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with the degeneracy factor g ~ 2 (left-handed neutrinos and right-handed an- 
tineutrinos). Equs.(0.1l and (0.2 I lead to 
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Introducing the (dimensionless) variables p and x through 
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equ.(0.3) becomes 
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The variable p is related to the gravitational potential V{r) through 
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Vq is the gravitational potemtial at the outer part of the galaxy, where the 
density vanishes. Indeed the numerical solution of equ.( |0.8| always provides 
a finite xq, where p{xq) — 0. The radius R of the galaxy is then R = qx^ 
(equ.(0.6)). The total mass of the galaxy (M = Mc + M^) is related to the 
derivative of p{x) at xq by 
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At r = Tc the boundary condition is 
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Equation (0.8 1 together with the boundary conditions (0.10 1 and ( 0.11[ ), repre- 
sents a loaded polytrope |H]-[II]. Similar mathematical structures appear in the 
study of a galactic nucleus, a neutron star, or in the Thomas-Fermi description 
of an atom. 
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We solved numerically equ.(0.8l, using the boundary condition of equ.(0.11l 
and an arbitrary positive value for the derivative of p{x) at x = Xc- The nu- 
merical solution provides a function p{x) which rises up to a maximum, then 
decreases until it vanishes at some point xq {p{xo) = 0). We evaluated also the 
derivative of p{x) in the neighbourhood of Xq and subsequently the left-hand 



side of equ.( |0.10 1, thus determining the ratio M^/Mc- To obtain the desired 



ratio Mt,/Mc, we repeat the numerical evaluation with a different value for the 
derivative of p{x) at Xc, until we achieve the predetermined M^/Mc ratio. The 
numerical study revealed the following essential features of the loaded polytrope: 

1. The precise value of xq is largely independent of the ratio M^/Mc- More 
massive neutrino halos provide larger values for the maximum oip{x), but 
all neutrino densities vanish in the vicinity of Xq — 2.0. Thus the radius 



R of the galaxy is set up by the constant q (equs.(0.6l and (0.7l) 



2. The scale q, equ.(0.7l, is very sensitive to the mass of neutrino m. By 
increasing the mass of the neutrino by a factor 3, the size of the galaxy 
is reduced by a factor 19. The overall data are well reproduced with a 
neutrino mass at 10 eV. 



Equ.(0.7l gives then the numerical expression 
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For a massive core Mc = 10^^ Mq, we obtain q ~ 3QKpc and therefore 
i? « 70 Kpc, in agreement with the data for the size of massive galaxies 
|12| . It is impressive that a scale q which is expressed in terms of funda- 
mental constants, such as the Planck constant, the neutrino mass and the 
Newton constant, reproduces accurately the galactic sizes. 



Through equ.(0.5l we obtain back the neutrino density N{r). We observe 



that at small distances, N{r) behaves as 
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While the neutrino density diverges as r -> 0, the total neutrino mass Af^ 
remains finite. 

Fig. 1 shows the neutrino density as a function of the rescaled distance x. The 
upper curve corresponds to Mc = 10^^ Mq, Af^ = lOAfc, while the lower curve 
corresponds to Mc = lO^^Af0, AI^, — Ale- Near the galactic core, the neutrino 
density is high as 10^ neutrinos/cm^, and over the uniform density of the big- 
bang cosmology, the gravitational clustering provides locally an increase by a 
factor 10^. 

The spherical neutrino halo up to a rescaled distance x, gives a mass fi^ix), 
where 

ti,{x)^aMQj yi/2[p(y)]3/2dy. (0.14) 

Obviously M^, = (i^^xq). The galactic rotational velocity, due to the neutrino 
halo, is then 



(0.15) 



3 



Fig. 2 shows the galactic rotation curves (upper curve corresponds to Mc — 
1O"M0, = lOMc, lower curve corresponds to Mc = lO^M©, = Mc) 
and there is agreement with the trends of the experimental data. 

We may reach our findings by making appeal to an even simpler model [l3]. 
A gravitational potential well of the size of the galactic halo is filled up with 
neutrinos, which are treated like a Fermi-Dirac gas at zero temperature. The 
total number of neutrinos is related to the Fermi momentum, itself determined 
by the minimum kinetic energy required for the neutrino to escape from the 
galaxy. It is found then that the radius of the galaxy is i? = g3^/'^, in full 
agreement with our results. 

Is it possible to detect the cosmological relic neutrinos? The gravitational 
clustering effect we studied here, enhances the local densities, however the aver- 
age neutrino energy is very small and the corresponding weak cross-section (of 
the order of lO""'^'^ crn^) renders detection by conventional means highly unlikely 

The few eV mass scale for the sterile neutrino, that our study suggests, 
fits nicely with the findings of the reactor and short-based neutrino oscillation 
experiments [5]. On the other hand the cosmological verdict for a few eV-mass 
sterile neutrino is rather unclear. The standard ACDM framework provides the 
constraint m < 0.5 eV. A modified though ACDM cosmology, with additional 
radiation, can accommodate a few eV neutrino|15j. In another direction, there 
is strong evidence for the existence of a 17 keV massive neutrino[TB]-|T7]. The 
coexistence of distinct mass scales for the neutrinos, sub-eV, eV and keV, indi- 
cates the complexity of neutrino physics. Clearly we need further experimental 
information and novel theoretical insights to decode the hidden dynamics. 

In summary we presented a simplified model for the gravitational clustering 
of massive neutrinos in the galaxies. There is no adjustable parameter in our 
analysis and despite its simplicity the overall features of galactic dynamics are 
reproduced. 
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Figures Caption 

Figure 1 : The neutrino density in a spherical galaxy as a function of the rescaled 
distance x. Upper curve corresponds to M^, — lOMc, the lower curve cor- 
responds to All, = Mc- For a Mc = IO^^Mq, real distance is obtained by 
multiplying the dimensionless x by 36 Kpc. 

Figure 2: The galactic rotational velocity, due to a spherical neutrino halo, as 
a function of the rescaled distance x. Upper curve corresponds to M^, = lOAfc, 
lower to Mjy = Mc- For a Mc = lO^^M©, real distance is obtained by multiplying 
the dimensionless x by 36 Kpc. 
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